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Abstract 

By remaining strictly within the confines of QCD, i.e., without invoking the 
parton model or any other exogenous element, we identify and resum pertur- 
batively to leading order, along with a correction term, enhanced radiative 
gluon contributions to the Drell-Yan type (qq pair annihilation) and deep- 
$h ' inelastic-scattering type (eq — > eq + X) cross-sections. The key feature of 

a ; 

the adopted approach is the recasting of QCD in terms of a space-time mode 

of description, which employs a path-integral formulation of field theories, as 

originally implied in works of Fock, Feynman and Schwinger. 
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I. INTRODUCTION 



Studies surrounding the infrared (IR) behavior of transition amplitudes in gauge field 
theories present special interest not only for the fundamental - of interpretational nature - 
questions they pose but also for the calculable physical implications they entail. To become 
more specific and, at the same time, identify the object of interest in the present investi- 
gation, we are here alluding to resummation procedures of perturbatively enhanced pure 
QCD contributions to physical cross-sections. These matters have come under systematic 
scrutiny by Sterman and collaborators |T|-[3|], following earlier realizations M on the physical 
relevance of soft gluon radiation. 

The wide acceptance of QCD as the microscopic theory of the strong interactions, despite 
the fact that its fundamental field quanta do not register as observed asymptotic states, 
rests largely on one's ability to isolate and identify its contribution to a given hadronic 
process by means of factorization theorems. With respect to the impact of IR properties 
of QCD on measurable quantities, the studies conducted in 0-0, which commence at the 
hadronic level and methodically proceed to the microscopic one |J, encourage the idea 
to pursue IR issues through a direct reference to QCD, i.e., through considerations which 
base themselves on the Green's functions of the theory. Recognizing, at the same time, 
the key role played by the parton model in connection with applications of QCD, we shall 
employ for its description a Fock-Feynman-Schwinger or equivalently worldline casting 
H HH (f° r recent reviews, we refer to [|l2i[l3| ), wherein a space-time interpretation - as 
opposed to a Feynman diagrammatic one - of fundamental matter field quanta propagation 
is attained flifl . In this way, one aspires to achieve consistency with the parton model 
mode of description without, however, invoking its probabilistic content. Our efforts, in this 
paper, will be directed towards determining pure QCD contributions related to processes 
induced by electroweak interactions - Drell-Yan (DY) and deep inelasting scattering (DIS), 
in particular. 

For the IR domain of QCD, a central issue of concern is where the perturbative picture 
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ends and to what extent non-perturbative effects start to become dominant (for a review 
and references, we refer to ||15|| ). As a matter of fact, the very concept of soft gluon radiation 
looses its meaning below some momentum scale of the order of Aqcd (recent discussions of 



these issues are given in [I6|,r7|). One consistent approach for attacking this problem strictly 
within the framework of the microscopic theory itself has been articulated by Ciafaloni |TB[ . 
As already mentioned, our own considerations will rely upon the worldline casting of QCD, 
keeping a more pragmatic course and focusing our attention on cross-section expressions. 

One might reasonably wonder as what constitutes the innovative tool, provided by the 
worldline casting of QCD, that facilitates calculations which "flirt" with the edge of the 
perturbative domain of the theory. A direct answer can be given at this point though its full 
content will become more transparent through the main exposition in the sections to follow. 
In short, it is the ability to isolate a special set of space-time paths having a very simple 
geometrical profile shared in a restricted (but directly relevant to the physics of the process) 
neighborhood by each and every contour entering the path integral. The single (multiplica- 
tive) renormalization constant carried by this special family of paths automatically factorizes 
their contribution to amplitudes/cross-sections given that it also accompanies the rest of the 
paths. The more complex geometrical structure of the latter simply implicates additional 
ultraviolet (UV) singularities which can be absorbed into conventional wave- function and 
coupling-constant renormalizations. This clean, geometrically based, argument, which sin- 
gularly underlines the worldline description, not only constitutes a notable simplification 
over the type of reasoning that has been promoted by Collins et al. || within the Feynman 
diagram perturbative approach to realize factorization in QCD, but, potentially, it paves 
the way to further applications. We especially have in mind issues pertaining to projections 
of our approach into the non-perturbative domain of the theory. 

Letting these comments suffice for an introductory exposition, we now proceed to display 
the organization of this paper, which is as follows. In the next section, we exhibit the world- 
line expression for the full fermionic Green's function and subsequently employ it to construct 
corresponding expressions both for DY and DIS type QCD amplitudes/cross-sections. Sec- 



tion 3 furnishes, with the aid of an Appendix, our basic calculations associated with one 
virtual gluon exchanges for a special set of trajectories. The resulting expression explicitly 
reveals the threshold enhancement factor while the task of virtual gluon re- summation is 
performed, via the aid of the renormalization group, in Section 4. Section 5 deals with the 
re-summation of contributions from real gluons, whereas our conclusions are presented in 
Section 6. 



II. BASIC WORLDLINE EXPRESSIONS FOR AMPLITUDES AND 

CROSS-SECTIONS 



Consider the full two-point (fermionic) Green's function in the presence of an external 
gluonic field. The expression, in Euclidean space-time, 



iGij(x,y\A) 



dTe~ Tm2 /«(„)=« Vx(t) 



x exp 



4 Jo 



x(T)=y 

dtx 2 {t) 



m — -7 • x(T) 



Pexp 



dtx ■ A(x(t)) 



displays the basic worldline features pertaining, more generally, to n-point Green's functions^ 
and, by extension, to amplitudes. Here, and below, P denotes the usual path ordering of 
the integrals. The first thing to point out is that a given path of the matter field quantum, 
starting at x and ending at y between respective "proper-time" values and T, also enters 
a Wilson line factor. The latter, being the sole carrier of the dynamics, separates itself from 
the rest of the factors in the path integral which are associated with geometrical properties 
of paths traversed by spin-1/2 particle entities. The most notable such quantity is the so- 
called spin factor |19[, Pexp (i/4) J T dta ■ 00 , where 00^ = (T/2)(x II x u — x fl x u ), accounting, 
in a geometrical way, for the spin-1/2 nature of the propagating particle. Accordingly, our 
perturbative expansions should be perceived of in terms of (Euclidean) space-time paths 



In this paper we shall be exclusively concerned with electroweak vertex functions of quarks which 
include gluon radiative effects. 
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involving a "proper time" parameter and not in terms of Feynman diagrams. As it turns 
out ||20|| , in the perturbative context, the structure of matter particle contours, entering 
the path integral, is determined by the points, where a momentum change takes place, i.e., 
points where a gauge field line (real or virtual) attaches itself on the (fermionic) matter field 
path. The almost everywhere non-differentiability of these contours, is residing precisely at 
these points. A major effort, in this paper, will be devoted to the extension of the worldline 
formalism to expressions for cross-sections corresponding to the particular processes of qq 
annihilation and e + q — > e + q + X . 

From the worldline point of view, the processes we intend to study involve fermionic 
matter particle (quarks) paths that commence at x and end at y, being forced to pass 
through an intermediate point z, where a momentum transfer Q takes place. This means 
that the Green's (vertex-type) function we shall be dealing with has the following form (r M 
denotes some Clifford-Dirac algebra element) 



V^ij(y,z,x\A) = GutfazlAlTpGkjfaxlA) 



where 



f dTe Tm2 [ ds [ x(0 )=x T>x(t)5 (x(s) - z) V M (x, s) exp 

JO JO J x (T)=v 

xP exp 



>x(T)=y 

r T 

ig / dixit) ■ Aixit)) 
Jo 



1 f T 



4 Jo 



dtx 2 (t) 



(2) 



J 'j 



r lt (x,s) = 



m — -7 • x(T) 



Pexpl-y dt a ■ uj I P exp I- J dta-uo 



m — —j ■ x{Sj 



(3) 



It is especially important to realize that in our approach off-shellness is naturally parame- 
terized in terms of the finite size of the matter particle contours and realistically accounts 
for the fact that quarks reside inside a hadron (m can be viewed as an effective quark mass). 
Going over to momentum space, we write 



V^j{p,if\z\A) = / dTe 
J 



-Tm 2 



T 



ds / Vx(t)S (x(s) — z) r M (i, s) 
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l r T 



4 Jo 



x exp 
x P exp 



= ^f[C 2 ]Pexp 

c z 



dtx 2 (t) + • x(0) + ip' ■ x(T) 



ig / dtxif) ■ A(x(t)) 



ig dx ■ A(x) 



(4) 



where C z denotes a generic path forced to pass through point z, at which the momentum Q 
is imparted. 

For a process of the type (DY) q + q — > lepton pair + X the "amplitude" expression 
reads 



.(DY) 
AW 



s')(— «7 ■ p' + m)V fi> ij(i , y ■ p + m)u(p, s) 
E^^l^expL/ dx-A(x) 

qz L JC Z 



(5) 



with the second, comprehensive, expression to be understood having recourse to Eq. 
Similarly, for the DIS-type process e + q — > e + q + X, one writes 



A 



(DIS) 
AW 



u(p', s')(«7 ■ p' + m)V^ij{i^ ■ p + m)u(p, s) 



Z(7 / <ix • 

Jc z 



(6) 



with the replacement p' — > — p' to be also made in Eq. (|J). 

For the cross-section, we need to employ the following quantity, which we implicitly 
display in Minkowski space-time after straightforward adjustments, 



C z ' c z 

xTr <P exp 



ig / dx" i A ll (x') 
Jc z 



Pexp 



(7) 



where P denotes anti- path ordering and the index a stands for either DY or DIS. Even 
though not explicitly displayed, the cross-section acquires a path-integral form, which has 
the following characteristics: 1) Paths C z and C z ' are forced to pass through points z and 
z' , respectively, where the momentum transfer occurs (see Fig. [I]). The distance b = \z — z'\ 
serves as a measure of how far apart the two conjugate contours can venture away from each 
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FIG. 1. Illustration of two conjugate contours C and C entering the worldline path integral, "talking" 
to each other at points z and z', where the momentum transfer for the physical process takes place. The 
distance \z — z'\ is referred to as the impact parameter. 

other and will be referred to as the impact parameter. 2) The traversal of C z ' is made in the 
opposite sense relative to C z , pretending the two paths join at one end (using translational 
invariance) , while allowing the other two ends of the contour to close at infinity, giving rise to 
the formation of a Wilson loop. 3) Under the circumstance just described, the Wilson loop 
formation guarantees the gauge invariance of the expression for the cross section. On the 
other hand, by keeping the contour lengths finite, but very large, thereby placing the quarks 
off-mass-shell, gauge invariance will still continue to hold to the order of approximation we 
employ in our computations, given that the off-mass-shellness also serves as an IR cutoff. 

Up to this point our considerations have been centered around the geometrical profile 
of the worldline casting of QCD. In terms of our comprehensive formulas @ and [FD for 
the amplitude and cross-section, respectively, the primary emphasis has been placed on the 
first factors of the respective sums, the main conclusion being that the relevant contours 
entering the path integral are open for the amplitude and closed (or almost so) for the cross- 
section. Armed with this information, we now turn our attention to the Wilson factor which 
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contains all the dynamics of the given process. The obvious task in front of us is to assess 
its implications once the gauge fields are quantized, i.e., once the Wilson factor is inserted 
into a functional integral weighted by the exponential of the Yang- Mills action. We display 
the quantity of interest as follows 



W = (TriPexp 



ig I dx'^AJx' 

< C z< 



A 



Pexp 



-ig J dx^A^(x) 



= (Tr(U\C z ')U(C z ))). (8) 

In the above expression, {■ ■ -}a signifies the expectation value with respect to the gauge 
field functional integral which, in this work, will be considered in the context of perturbation 
theory. Note in the same context that a virtual gluon attaching itself with both ends to the 
fermionic worldline, entering the amplitude, corresponds to a correlator between a pair of 
gauge fields originating from the expansion of the Wilson factor. On the other hand, for an 
emitted "real" gluon from the fermionic line, the correlator is between an "external" and 
a Wilson-line gauge field.Q The overall situation is depicted in Fig. |^. At the cross-section 
level, now, "real" gluons are integrated with respect to "propagators" linking together the 
two conjugate contours, while their polarization vectors are summed over (cut propagators). 
This is precisely what (• • •) signifies in the last equation, as it covers both Wilson line factors. 

In the light of the above remarks, let us proceed to display the first-order (in perturbation 
theory) expression for W, which receives contributions both from virtual gluons, viz.; ones 
attached at both ends either to worldline contour C z or to contour C z> , as well as from 
"real" gluons linking these contours to each other (cf. Fig. [|). This expression reads 

W (2) = Tr / - g 2 C F F dt x f dt 2 6 (t 2 - t x ) a? (t 2 ) x v {t x ) (x(t 2 ) - zfa)) 
Jo Jo 

- g 2 C F [ T dt[ [ T di! 2 9 (ti - t' 2 ) x ,fl (t' 2 ) x w (ti) Dau {x(t' 2 ) - x(ti)) 

Jo Jo 

- g 2 C F £ dt £ dt' x(t) ■ x' (f) D cnt (x(t) - x(t')) + O (g 4 ) . (9) 



2 One will, of course, also encounter correlators that involve gauge fields from the non-linear terms 
of the Yang-Mills action. These, however, do not enter the leading logarithmic considerations. 
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c 



x(s') = z 




C 



z = x(s) 



x(0) ' 1 1 x(0) 

FIG. 2. Virtual gluon radiative corrections of various sorts and "real" gluon lines with their ends 
attached on each of the two depicted contours at the cross-section level. 

It becomes obvious from their structure that the first two non-trivial terms correspond to 
virtual gluon contributions - one per conjugate branch -, while the third one is associated 
with "real" gluon emission. Finally, concerning the gluon propagators entering the above 
equation, we shall be employing their Feynman-gauge form without loss of generality due 
to gauge invariance. In particular we have, in .D-dimensions, 



whereas 



A-D 



d D k 



-ik-x 



(2tt) d k 2 + tO. 



9fj.u 



4tt 2 



-7T/i' 



(4-z>)/2 T(D/2 - 1) 



( x 2 _ iQ+ )(0/2)-l 



(10) 



A-D 



d u q 



D 



27rd(q 2 )6(q°)e 



0\ — iq-x 



4tt 2 



-7T/i 



a \(4-Z>)/2 



V{D/2-l) 



[{xl - *0_ 



2 ,(2J/2)-l- 



(11) 



As already established by other methods, the perturbative expansion (||) is plagued by 
large threshold logarithms leading to the need for factorization and resummation. This is 
precisely the task we are about to undertake within the framework we have developed and 
which does not take us outside QCD. 
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III. FIRST-ORDER VIRTUAL GLUON CORRECTIONS UNIVERSAL TO THE 

VICINITIES OF POINTS z AND z' 



The space-time mode of description of our worldline casting puts us into the position to 
promote the following argument: The point z (or z'), where the momentum transfer Q is 
imparted, marks the presence of a neighborhood around it, no matter how infinitesimal in 
size this might be, whose geometrical structure is shared by all fermionic paths entering the 
path- integral. Specifically, there will be a derailment (cusp formation), whose opening angle 
will be fixed unambiguously, since it is determined by the momentum transfer. It follows 
that the contributions to the amplitude and cross-section from the immediate vicinities of 
these points is a common feature of all contours and eventually factorizes. In this section 
we shall determine the first-order perturbative term corresponding exactly to this factor. 

Consider now the neighborhood of point z on the contour C z . Expanding around this 
point, we write 

x»(t) = af(s)+tx» l (s±a) + ... (12) 

with v M = x M (s — 0) and v ' M = x M (s + 0) being entrance and exit four- velocities, respectively, 
with respect to z. 

Adjusting our notation by re-parameterizing the contour so that the zero value is assigned 
to point z, the relevant quantity to compute, to first perturbative order, becomes 



Ugl = l- 9 2 Q 



f° dt x [° dt 2 9 (t 2 - ti) u'VD™ {vt 2 -vti)+ r dt x j" dt 2 9 {t 2 - t x 

J-a J -a ' JO JO 



+ vVD^ (v't 2 - v't x ) + J " dh £ dt 2 6 {t 2 - h) v'WD^ {v% - vt x ) 



(13) 



It is clear that the above expression corresponds to the first term in Eq. (|9]), which monitors 
virtual gluon exchanges occurring on contour C. 

From the above considerations it follows that the main contribution to each double in- 
tegral comes from the common limit tx, t 2 — ► 0. Suppose now, the other limit is to be 
determined by demanding that its contribution to the integrals is of vanishing importance. 
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Then, such a requirement automatically isolates those contours, whose only significant ge- 
ometrical characteristic is that the four-velocities to approach and depart from point z are 
fixed, denoted by t> M and v ,fl respectively; the same, of course, happens for point z', but 
in reverse order. This justifies the subscript S in Uc,Si which stands for "smooth" . More 
complex geometrical properties such as non-differentiability almost everywhere (which most 
certainly characterizes the vast majority of paths) do not enter Eq. ([13]). The point is, on the 
other hand, that every path will share the geometrical structure we are focusing on in some 
neighborhood of the point z, or z', no matter how close to these points one has to come. 
At the same time, the UV singularities of this restricted set of paths will entail expressions 
that solely depend on the two four- velocities and the opening angle. Paths of more complex 
geometrical structure, on the other hand, will certainly exhibit these UV singularities plus 
additional ones.0 It naturally follows that our restricted set of trajectories, by exclusively 
determining the corresponding (multiplicative) renormalization constant will factorize from 
the rest of the expression for the amplitude and/or the cross-section. From dimensional 
considerations it follows that the omitted terms in Eq. QT2D will contain negative powers of 
a, whose dimension in the denominator is (mass) 2 . Neglecting their presence means that 
a should be very large in magnitude and hence should be related to the IR cutoff, i.e., 
a ~ A~ 2 , A > Aqcd- 

Going over to Minkowski space-time, there are two distinct possibilities for defining an 
infinitesimally small neighborhood around z. The first one, to be labelled (a), reads 

(x-x'f = C(e 2 ), with ~ uj,, for all //, (14) 

where e(< Q^ 1 ) is a small length scale. The second alternative, to be labelled (b), can be 
typically represented by 



3 Actually, the standard UV singularities of perturbative field theories associated with j3- 
functions, coupling constant and wave-function renormalization pertain to almost everywhere non- 
differentiable paths. 
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(x - x') 2 =0(e 2 ) with \v - v'\ 2 = 0(X 2 ) but (y+ - v' + ) ~ O(Q) and (v_ - ?/ ) ~ C 

=►(«+-<)(«_ -t/_) = 0(A 2 ) (15) 

that is equivalently effected via the condition v + ^> v' + , V- ~ v'_. All in all, there are four 
different configurations: + <-» — and prime^no-prime entering this case. 

We characterize case (a) as "uniformly soft" , given that the considered gluon exchanges 
take place in a neighborhood whose smallness pertains to all directions. Case (b), on the 
other hand, will be referred to as "jet" since gluon emission occurs under circumstances, 
where entrance and exit four-velocities differ between them significantly along one or the 
other of the light-cone directions. We point out that what here is termed as "jet" does not 
coincide with the notion of Collins et al. ||. The actual connection will become evident 
later on. 

Let us commence our calculations by taking up the first 0(g 2 ) term entering the right 
hand side of Eq. ([TBI). Since this only involves the branch of the contour C z entering point 
z, we obtain the same expression regardless of whether or not a uniformly soft or a jet 
configuration is being considered. It reads 



h = [ dti f dt 2 6 (t 2 - h) v»v v D^ (vt 2 - vh) 

J —a J —a 



1 '-^r^-^^. <h» 



8tt 2 V ^ v V 2 / D-3 2- D/2'' 

where^ L\ = u\v\. The second term has the same structure (it involves the exiting branch 
of C z ) and therefore produces a similar result: 

2 8tt 2 V p V \2 J D-3 2- D/2 K J 

with L 2 = cr\v'\. 

A couple of remarks are in order at this point. First, even though the length scales L\ 
and L 2 are both large, being proportional to a, they will be of the same order of magnitude 



4 Note that v has dimensions of mass as our "time" parameter a has dimensions of (mass)" 
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for case (a), whereas for case (b), one scale will be negligible in comparison with the other. 
Accordingly, the total expression for the uniformly soft amplitude will be twice as large 
as that of the jet-like one. This being said, we shall denote the dominant length scale by 
L(~ L\ and/or L 2 ), when it enters our final expressions, and set it equal to j, recognizing 
that it is of the same order as the IR cutoff. Second, in order to avoid the double counting 
resulting from the fact that each branch has been "cut-off" at corresponding distances away 
from z (implying that gluon emission will occur at the endpoints to be offset by a similar 
one but opposite in sign from that portion of the contour that continues to stretch out to 
infinity), the final expressions for the amplitudes/cross-sections should be multiplied by a 
factor of 1/2. Equivalently, one might think of this compensation as actually identifying 
the missing energy of the gluon emission at the extremities with the off-mass-shellness. In 
fact, this is what we have been implying all along when claiming that finite contours signify 
off-mass-shellness. 

Turning our attention to the contribution resulting from a virtual gluon exchange from 
the entrance to the exit branch, with respect to z, we consider the quantity 

/0 ra 1 , \(4—D)/2 /D \ 

^ dt, j o dt 2 V'^D^ ( V 't 2 - Vh) = — (-7T/! 2 ) » r(--l) 

x v ■ v' J dtx J dt 2 (t\v 2 + t\v' 2 + 2t x t 2 v ■ v' - i0 + J . (18) 

For case (a) it assumes the form (recall that v • v' is negative for the DY and positive for the 
DIS type of process) 

x[ dt 2 lt 2 + t 2 2 + 2t 1 t 2 ^--i0 + ) . (19) 

Jo \ \ v \\ v \ / 

As shown in the Appendix, for DY, one determines (7e is Euler's constant) 

j Sy = ^^ DY coth ^ y ^d; + ^^ DY coth ^ DY ln (^ 7re2+7E ) ' ^ 20 ) 

where cosli7 DY = w DY = > 1, while for DIS (w DIS = > 1) 

13 



4%is = ^~2^ DIS coth7 D is 2 -j g + ^7dis coth7 D is In ^7re 2+7E ^J . (21) 

In all of the above expressions, as well as in those that will follow, we have ignored: (i) 
all imaginary terms that will drop out when contributions (for virtual gluons) from the 
conjugate contour are taken into account and (ii) finite, /i-independent terms that will 
cancel out when real-gluon contributions to the cross-section are included. 
Collecting all terms, we deduce, for the "uniformly smooth" part, 

li a) + 4 a) + 4 a) = ^(7coth 7 - 1)^ + ^( 7 coth 7 - 1) In (^e 2+ ^ , (22) 

where the angle 7 is to be adjusted either to the DY or the DIS situation. 

Concerning the "jet" part of the computation, we only need to consider I3 because | 
the expression for 7^ + 7^ is simply one half of that of l[ a ^ + ■ A typical term entering 
4 b) 0- > V+) is 

1 / 9 n(4-d)/2^ (d \ , r ,, r /.9 9 . „. . , .„ \i-(^/2) 



J<» = -±- 2 (-vr/i 2 ) (4 ^ )/2 r (I - l) „ ■ v> £ dt, £ dt 2 {t\v 2 + 2 tl t 2 v -v'-iO 



(23) 



whose computation suffices to furnish each of the other three terms as well. 

It is shown in the latter part of the Appendix that for either the DY or the DIS type 
process one obtains 

& = is? fZff + «b ?V (f- ,e ) + sh ln2 (iH + < 24 > 

Subtracting the pole terms in the MS scheme, we arrive at the finite part of the overall 
result. For the uniformly soft contribution, in particular, we get 

(7{ a) + + 7^) fin = JL ( 7Cot h7 - 1) In (tpj , (25) 

while the jet contribution reads 



5 Recall the remark following Eq. (|l 
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(4» + 4* + 4V = ^>(£), (26) 



where we have set A 2 = 4A 2 e _27E . The above relation takes into account all four different 
configurations contributing to I^ b \ 

Putting everything together we arrive at the following overall result for the second order 
contribution stemming from contour C z 



7 coth 7lnf^Tr| lnf^rr) +ln 2 f^TT 

UV 2 UV Va 2 , 



(27) 



2tt 

A similar result is obtained also for contour C z . 

Noting that 7coth7 = ln(Q 2 /m 2 ) (for Q 2 3> m 2 ), with Q 2 = (p + p') 2 for the DY 
and Q 2 = —{p'~p) 2 for the DIS type of process, respectively, we recognize that the well 
known perturbative enhancements occurring as Q 2 — > oo are associated with the eikonal-type 
trajectories upon which our present calculations have been based. One, now, realizes that 
these trajectories define threshold conditions, with respect to the given momentum exchange 
Q, for the processes under consideration since they leave no room for space-time contour 
fluctuations. In the following section we shall deal with the summation of these enhanced 
contributions to leading logarithmic order. We shall, furthermore, identify a correction factor 
associated with those terms in Eq. fl27|) not involving the enhancement factor ln(Q 2 /m 2 ). 



IV. SUMMATION OF ENHANCED CONTRIBUTIONS FROM VIRTUAL 

GLUONS 

The family of worldline paths to which the considerations in the previous section have 
referred to recognize all (virtual) single-gluon exchanges, consistent with the simple geomet- 
rical configuration of two constant four-velocities making a fixed angle 7 between them (in 
Euclidean formulation). Among these gluons there will be "hard" ones (upper limit Q) and 
"soft" ones (lower limit set by A). What is debited to the former and what to the latter 
group of gluons is, of course, relative. It is precisely the role of the renormalization scale 
/i, entering through the need to face UV divergences arising even for the restricted family 
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of paths, to provide the dividing line. The corresponding renormalizat ion-group equation 
reflects the fact that the scale \x is arbitrary and that physical results do not depend on it. 
A straightforward application of this fact will enable us to sum the enhanced, virtual gluon 
contribution to the amplitude in leading logarithmic order as well as to obtain a bona-fide 
correction term. 

To bring the above discussion into a concrete form, consider a factorization, good to 
order 1/Q 2 , between the "soft" and the "hard", relative to the scale /i, gluon contributions 
to the amplitude Uc entering Eq. (U) for the restricted set of contours.0 We write 



Uc = U c A$4)ucM) + o(^) (28) 



(2) 

noting that the soft part Uc,s coincides, to second order, with Uq S . 

It is further convenient to separate the part of the "soft" term involving the enhancement 
factor from the one which does not. To this end we write, to second order in perturbation 
theory, 

Uc,S = ^C.cusp^C.coll (29) 



with 



^L p =l-^.n(S)ln(£) <30) 



and 



U C,coll — 1 _ 27T F 



2\ Q / ,X 

(31) 



k2 ? -H? 



In the above relations, the designation "cusp" refers to that part of the soft factor which 
recognizes the angle 7, whereas the characterization "collinear" signifies independence from 
7. From a physical standpoint, it represents the contribution from those gluons which 



6 For the purposes of the following discussion, it suffices to use the generic notation C for either 
contour. 
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are emitted collinearly with incoming or outgoing four-velocities due to the sudden accelera- 
tion/deceleration occurring at the cusp point. It should be further pointed out that collinear 
emission occurring in the above context pertains to the soft sector and does not coincide in 
the full with what is designated as "collinear" by Collins et al. in ||. Full agreement on 
"collinearity" will arise a posteriori, i.e., once resummation is performed. 

From Eqs. ([28l)-(l31~|) we obtain a relation which does not involve U co \\ by simply taking 
the logarithmic derivative with respect to Q 2 : 



d 



In U = -rA^ In U cusp + - In U u , 



(32) 



d\nQ 2 d\nQ 2 LUbp d\nQ 2 

where we have dropped the subscript C as being superfluous in the considerations to follow 
The /i-independence of physical results leads to the renormalization group equation 



d 



d 



cHn/i d\nQ- 



ln% = -^rf!nV lnf/c - 



2 1 cusp l 



with r cusp to be read off from Eqs. (E2])-(E|) and (O); specifically, 



(33) 



r C us P (a s ) = -c F + oK). 

7T 



From the second leg of Eq. fl3"3|), one obtains 



dlnQ 



2 In t^c US p 



v 2 dt 

A* 2t 



r c us P [a 8 {t)] 



which, in turn, gives 



d 



d\nQ 2 
where we have defined 



In Ok 



Q t p [a s (i)]+rL(Q 2 ) 
p? zt L 



(34) 



(35) 



(36) 



a s (Q 2 



rflnQ 2 

Combining the last three equations we have 



\nU H 



V / p 2 =Q 2 



d 



d\nQ' 



In £7 



Q 2 fa 

_ -r cusp [a 8 (t)} + r a s (Q 2 ) 

X 2 Zt 



(37) 



(38) 
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where, in terms of In U , we write 



« S (Q 2 



dlnQ- 



\nU\ 



(39) 



Setting Q — fj, , we are led to the isolation of the "cusp" and "collinear" terms according 



to 



,„ U = in U„ in (f ) + .„ ^ ^ 

from which, once appealing to Eq. (|35D , we obtain 



d 



d\nQ 2 

which leads to the identification 



\nU\ 



X 2 =Q 2 



d^ lnUc0lll ~ X2 = Q2 



(40) 



(41) 



« S (Q 2 



In C/coiil A2= Q 2 = -MQ 2 ) + 0{a 2 s ). 



d\nQ 2 ""^ LL,A ~ =w ~ 4 



(42) 



Collecting together all our findings, we obtain our final, re-summed result corresponding 
to the contour C. It reads 



U 



exp 



Q 2 dt 
A2 2t 



■ Q2 

In— r cusp (« s (t)) - T(a s (t)) 



U (a s (Q 2 )). 



(43) 



One notes that the correction factor is associated with collinear emission (cf. Eq. (f42D ). 
Secondly, the factor Uo(a s (Q 2 )) represents initial conditions input at the QCD level. One 
cannot help but bring to mind evolution equations at the partonic level where the proba- 
bilistic interpretation prevails. The connection between evolution and renormalization group 
equations is widely taken for granted. The worldline approach might provide the space-time 
description framework for understanding the precise connection between the two types of 
equations. Finally, let us remark that the conjugate-contour term W(C Z ') can be treated in 
a completely analogous fashion. 
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x(s') = z' 



z = x(s) 



FIG. 3. Neighborhoods of respective points on two conjugate contours, where the momentum transfer 
takes place, and associated four-velocities. 

V. SUMMATION OF ENHANCED CONTRIBUTIONS FROM REAL GLUON 

EMISSION 



We shall now turn our attention to real gluons and attempt to factorize cross-section 
contributions from neighborhoods around points z and z' . Note that we now have to deal 
with gluons which connect two "opposite" neighborhoods while crossing the unitarity line, 
(this situation is depicted in Fig. The relevant scale promptly entering our considerations 
is the impact parameter b = z — z', which must be eventually integrated over in order to get 
the cross-section. Naturally, the short-distance cutoff in this integration will be provided by 
the (length) scale l/\Q\. 

For the eikonal-type family of paths and in first order perturbation theory, the relevant 
quantity on which our quantitative considerations are to be based, i.e., the counterpart of 



Eq. (13), is given by 



U, 



(2) 



cc,s 



l + 9 2 C F 



rO 

dt\ j dt 2 v ■ v D (tiv — t 2 v + b) 



I dt x \ dt 2 v' ■ v' D {t x v' - t 2 v' + b) 
10 Jo 
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+ f° dtx f 7 dt 2 v ■ v' D (txv - t 2 v' + b) 

J -a JO 

+ f dt x [ dt 2 v' ■ v D (hv' - t 2 v + b) 

JO J -a 



(44) 



where the bar denotes four-velocities for the conjugate contour and the subscript cut is 
henceforth omitted. 

(2) 

To identify the leading behavior of [/^ s , with respect to b, we shall consider first the 
situation corresponding to b = 0. The subsequent emergence of UV divergences, once han- 
dled through dimensional regularization, will introduce a mass scale p, that will be bounded 
from below by an IR cutoff A and from above by the (mass) scale 1/6. The resulting renor- 
malization group equation will facilitate the resummation of the leading terms, just as in 
the virtual-gluon case. 

Let us commence our quantitative considerations by looking at the term 

rO rO 

J 1 (6)=w-w/ dtx / dt 2 D (hv - t 2 v + b) (45) 

with v 2 = v 2 = —v ■ v (see Fig. |3|). 

Setting 6 = and using the expression for the cut propagator as given by Eq. (|TTD, we 
obtain 

m - -i 2 H**r» t (j- i )d^3t^d[ i - - «] <«) 

which actually coincides with what one would obtain if the regular propagator was substi- 
tuted. The significance of this occurrence is that it leads to the same anomalous dimensions 
for the running of the real gluon contribution to the cross-section as for the virtual part. 
This fact can be immediately verified via a direct comparison with Eq. (|T6|). 
Isolating the finite part of the above expression, we write 

■^ = -5? '"(§)' (47) 

It is trivial to see that the same result holds also for J 2 a ^'^ n . 
We next turn our attention to the term 
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J 3 (b) = v ■ v' I dt x r dt 2 D (ttv - t 2 v' + b) . (48) 

J-a JO 

Its computation will concurrently allow us to determine J^b) which corresponds to the 
exchange primes no-prime in the expression above. 

Dimensionally regularizing the cut propagator, we obtain 

£dt 2 {t\v 2 + t\v n + 2v ■ v't x t 2 - i0+) 1_D/2 . (49) 



x 



Once again we record, by referring to Eq. (p23|), coincidence of the singularities and, by 
extension, of associated anomalous dimensions between virtual and real gluon expressions 
that contribute to the cross section. 

For the "uniformly soft" configuration the corresponding result is 

l-D/2 



An 2 V 2 / If I If ' 



x [ dhf dt 2 lt 2 1 + t 2 2 + 2t 1 t 2 ^-^- -i0 + ) . (50) 
Jo Jo \ MM J 

Taking into consideration that for the DY case rjr=n = -tttft = cosh 7dy > we obtain 

° \ v \\ v \ v \\ v \ 

4 a DY (o) = JS Y (o) = ^(-V) (4 ^ )/2 r (| - l) COsh 7DY 



X 



J dt x j dt 2 (t\ + t\ + 2t x t 2 cosh 7dy - z0+) 1 °' 2 , (51) 
whose finite part reads 

jg Y n (0) = jSf(O) = -^ 7 DYCoth 7DY In (~Q . (52) 



For the DIS case, where = cosh(7 DIS - wr) < 0, and y^p = cosh (701s < 0, we 
find 

45is n (°) = g^a <^ DIS " i7T ) coth 7Dis In (^j , (53) 

whereas 
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J $is n (0) = (7dis + itt) coth 7DIS In . (54) 

Given that the imaginary parts in the last two equations will eventually drop out, from 
now on we shall write, generically, 

4 a)M (0) = 4 a)M (0) = ^7coth 7 In (^j , (55) 

where 7 can be adjusted to the appropriate case. 

Turning our attention to the "jet" configuration, we can go directly to J^ b \o), since 
j[ b \0) + J^ify furnishes half the contribution of its uniformly soft counterpart, the reason 
being the same as the one given in the virtual gluon case. We thus have 

(b), x 1 / 9\(4-£0/2 (D \vv' f 1 , f 1 , f o v-v' \ 1_D / 2 

J 3 l, (0) = r(--i) ¥ / ( *, jj| *, ( t f + 2M 2 ^ - io + j 

(56) 

with a corresponding expression holding also for J^ b \o). 
For the finite parts one obtains 

Jf fin (0) + Jf ) ' fin (0) = i ^ln 2 ^, (57) 

which holds true for both the DY and the DIS case. 

Collecting our findings from the real-gluon analysis to the second-order level, we write 
for the finite contribution to the cross-section 



UcCsw — H — ~Cf 

7T 



7coth7ln (g)_3 ln (a +ln2 (£ 



~ U CC,cub P U CCI,coU ' ( 58 ) 

where (to the order we have been calculating) the results for the "cusp" and the "collinear" 
terms read 

tWusp = 1 + ^C Fl coth 7 In (tj^J (59) 

and 
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u. 



(2) 

CC.coll 



7T 



-MS MS 



(60) 



At the same time, the singularity structure of the full expression for the cross-section 
entails a multiplicative renormalization factor, which is common to all "Wilson loop" con- 
figurations entering its description, but which is the only one that pertains to the family 
of eikonal-type paths under consideration. The reasoning is, of course, identical to the 
one promoted for the virtual gluon case. Therefore, the corresponding contribution to the 
cross-section factorizes and the same resummation procedure can be employed as for the 
virtual-gluon case. As already observed, the anomalous dimension is in both cases the same. 
There are, however, the following notable differences. First, the upper limit for the mo- 
mentum of real-gluon emission is 1/b 2 instead of Q 2 . Second, there is a difference of sign, 
which becomes evident by comparing Eqs. fl30"D and fl3~T|) with Eqs. (p9|) and ( j60| ) and, fi- 
nally, no hard real-gluon emission enters our considerations - by definition. In this light, it 
is practically self-evident that the resummed expression for real-gluon emission becomes 

"i/& 2 dt 



U CG = exp J 



i Q2 v 
In —I 



cusp 



(a s (t)) - T (a s (t)) 



U, 



cc,o- 



(61) 



We can now bring together real and virtual gluon results by referring to our generic 
expression for the cross-section as given by Eq. (H). We write 



W 



Tr (U\C Z ')U(C Z 



UcUcUcc, 



(62) 



where we have first extracted the resumed expression corresponding to the curly bracket ex- 
pectation values from virtual gluon exchanges in the amplitude and its conjugate. Clearly, 
the last factor corresponds to the resummed expression resulting from the real-gluon "aver- 
aging" . 

Now, at the cross-section level, our threshold resummation of the virtual gluons reads 



U c Uc = exp | - j 



Q 2 dt 

A2 T 



' Q2 

In— r cusp (a s (t)) - r(a B (t)) 



Uc,oU C n- 



(63) 



Finally, by making the re-adjustment In (/i 2 /A 2 ) — > In /\ 2 j we obtain our final result, 
which reads 
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y^cc = exp { - 



Q 2 dt 

/6 2 t 



with r cusp and T given by Eqs. 



t 



and 



cusp 



Wo, (64) 
, respectively and where c = 4e~ 27E . an 



expression whose previous derivation, see, e.g., Ref. |2l| and references therein, has employed 
Wilson lines, as an exogenous element attached to quark operators 



VI. CONCLUDING REMARKS 

In this paper we have studied the threshold resummation behavior of DY and DIS type 
of processes, widely pursued by Sterman and others (see, for example, in |]||), by staying 
strictly within the framework of QCD. From a physical viewpoint, we based our consider- 
ations on soft-gluon radiation, restricting ourselves to an energy regime whose lower cutoff 
is high enough to justify an analysis in which reference to "gluons", as dynamical degrees 
of freedom, continues to make sense. In this context we have implicitly assumed the pre- 
confinement property, originally articulated in the first work of reference Q (see also [Q), 
according to which the non-perturbative dynamics responsible for confinement screens color 
up to the infrared scale A which sets the lower limit for the perturbative regime. 

Our analysis has been conducted via the adoption of the worldline casting of QCD whose 
emphasis on a space-time propagation content facilitates the isolation and eventually factor- 
ization of IR contributions to cross sections through eikonal-type paths. In this connection, 
let us remark that we have already extensively applied this approach to the calculation of 



various fermionic Green's functions p3| (see, also the Sudakov form factor [25|, as well 



as to derive gluon Reggeization in connection with the four-point fermionic Green's function 



for the latter case in the forward regime [26]. The present paper presents the first applica- 



tion of the worldline methodology to "cross-sections" in which real-gluon radiation effects 
must also be taken into account - and so far have received much less systematic treatment. 

The common renormalizat ion-group running of virtual and real - uniformly soft and jet- 
type - gluon contributions to the cross-section has led to a final result in which the IR-cutoff 
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momentum is provided by the inverse impact parameter which can be vested with concrete 
physical meaning related to the process itself. A number of possible physical applications, 
ranging from extensions to exclusive processes to intrusions into the non-perturbative regime 



of QCD, as per, for example, along the lines of Ref. come to mind which will occupy 
our attention in the immediate future. 
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APPENDIX A: 



Our task is to establish Eqs. fl20|) and (21) in the text. Performing the integration 
entering the right hand side of Eq. (0), one obtains 

Aa) 1 / 2 \{±-D)/2 (D \ 1 1 



f / L> L>- 1 „ \ 

x {w DY F \1, — - 1; — ^— ; 1 - w DY J 



+^ [2(1 - wdy)] f ( 1 ^- i; V ; r (A1) 



for DY kinematics and 



' 3 ,DIS = ^ (-7TA* ) r f - - 1 j — — 2 WDIS 



.2 



x |wdi S F (JL, — - 1; — ^— ; 1 - w ms 

-\ [2(1 + » DIS )] 2 - D ' 2 F (l, £ - Izpfi) > (A2) 



for the DIS situation. 

Setting D = 4, we obtain 



F (l, 1; 3/2; 1 - ^ Y ) = — 7DY (A3) 

v ' smh 7dy cosh 7dy 



and 



F(l,l;3/2;i±f2l)=^- i -^. (A4) 
V 2 / smh7DY smh7DY 

As the imaginary part in the above expression will cancel against its counterpart in the 

conjugate expression, it can be dropped as far as the cross-section is concerned. 

Denoting the expression inside the curly brackets on the rhs of Eq. ( [All ) by fo Y (w) and 

setting 

fE Y M = /rn + (4 - V gY(w ]:f YH (as) 

we realize that the second term on the rhs will lead to finite terms that depend solely on 
w and which will cancel from similar contributions of the same sort coming from the other 
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terms entering Eq. ([13]). Putting everything together, one finally arrives at Eq. (p0|). Clearly, 
the DIS situation is similarly confronted. 

To establish the result given by Eq. fl2H|) , we first note that Eq. ( p3j ) gives 

, 2 \ (4-r>)/2 ,x (4-D)/2 



r(6) 



1 



4tt 2 



-7T- 



A 2 



2 / (4-D) 2 I M 2 



x 



D D D 2v-v'\ ( 2w' 
F(--l,2--;3--; — I +1 + 



2-D/2 



2 2'" 2' M 2 / V Iwl 2 



'2vv'' 



2-D/2' 



(A6) 



Now, the distinction between the DY and DIS like processes resides in the sign of v ■ v' 
(positive and negative, respectively). On the other hand, what is relevant for the "jet" 
situation is that the absolute value of 2j^p- is of the order of unity. One, then, easily 



determines that in the limit D — > 4 Eq. 



is retrieved. 
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